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ABSTRACT 

We present high-resolution spectroscopy of the neutron-star/low-mass X-ray binaries 2S 0918—549 
and 4U 1543—624 with the High Energy Transmission Grating Spectrometer onboard the Chandra 
X-ray Observatory and the Reflection Grating Spectrometer onboard XMM-Newton. Previous low- 
resolution spectra of both sources showed a broad line-like feature at 0.7 keV that was originally 
attributed to unresolved line emission. We recently showed that this feature could also be due to 
excess neutral Ne absorption, and this is confirmed by the new high-resolution Chandra spectra. The 
Chandra spectra are each well fit by an absorbed power-law + blackbody model with a modified 
Ne/O number ratio of 0.52±0.12 for 2S 0918-549 and 1.5±0.3 for 4U 1543-624, compared to the 
interstellar-medium value of 0.18. The XMM spectrum of 2S 0918—549 is best fit by an absorbed 
power-law model with a Ne/O number ratio of 0.46±0.03, consistent with the Chandra result. On 
the other hand, the XMM spectrum of 4U 1543—624 is softer and less luminous than the Chandra 
spectrum and has a best-fit Ne/O number ratio of 0.54±0.03. The difference between the measured 
abundances and the expected interstellar ratio, as well as the variation of the column densities of O 
and Ne in 4U 1543—624, supports the suggestion that there is absorption local to these binaries. We 
propose that the variations in the O and Ne column densities of 4U 1543—624 are caused by changes 
in the ionization structure of the local absorbing material. It is important to understand the effect 
of ionization on the measured absorption columns before the abundance of the local material can 
be determined. This work supports our earlier suggestion that 2S 0918—549 and 4U 1543—624 are 
ultracompact binaries with Ne-rich companions. 

Subject headings: binaries: close — stars: neutron — stars: individual (2S 0918—549) — stars: 
individual (4U 1543—624) — X-rays: binaries 



1. INTRODUCTION 

Low mass X-ray binaries (LMXBs) consist of a neu- 
tron star (NS) or black hole (BH) in orbit with 
a <1 Mq companion. An intriguing sub-class of 
LMXBs are the ultracompact binaries that have or- 
bital periods less than 80 minutes. Hydrogen-rich 
companions cannot sust ain a LMXB system with such 
a short orbital period <| Paczyn ski fc Sienkiewiczl Il981t 
iRappaport. Joss, fe Webbinkl Il98^i ~ However, orbital 
periods <80 min are predict ed for hydrogen-deficient 
or de generate companions llJoss. Ayni. fc RappaporH 
119781: iNelson. Rappaport. fc Jossl 119861) . and this was 
confirmed with the orbital period measurements 
of the X-ray pulsar 4U 1626-67 (P or b=42 min; 
iMiddleditch et al.|ll98l, the X-ra y dipp er 4U 1915-05 
(P nrh =50 min: IWhite fc Swankl IT982 IWalter et all 
1982) , the X-ray bursters 4U 182 0-30 (P orb =ll min; 
IStella. Priedhorskv. fc Wbrtl IT987h and 4U 1850-087 
(-forb=21 min: iHomer et al . 1996),' and the detection of 
the white dw arf analogs, the AM CVn systems (e.g., 
lWarnerl H~995'). In addition, the three recently discov- 
ered millisecond X-ray pulsars XTE J1751-305, XTE 
J0929-314, and XTE J1807-294 were also found to be 



LMXBs are the remains of He white dwarfs (WDs) 
that have transferred a significant fraction of their 
mass to the NS. However, recent X-ray spectral ev- 
idenc e indicates that s o me companions may be Ne- 
rich jSchulz et alJ 120011: l.luett. Psaltis. fc Chakrabartvl 
120011) . which suggests the possibility of C-0 or O-Ne- 
Mg WD companions. The growing population of these 
exotic ultracompact systems, as well as the new evi- 
dence for Ne-rich companions, has interesting implica- 
tions for the formation and evolution of binary systems. 
Motivat ed by the observational evidence for Ne-rich 
donors, lYungelson. Nelemans. fc van den Heuvell {2002) 
explored formation scenarios for these systems. Here we 
discuss two other NS/LMXBs that may also be ultra- 
compact binaries. 

Both 2S 0918-549 (I = 275?9, b = -3?8) and 
4U 1543-624 (I = 321?8, b = -6?3) have been ob- 
served by all of the major X-ray satellites since Uhuru. 
McClintock et al. (1978) identified the optical coun- 
terpart of 4U 1543-624 based on the SAS-3 posi- 
tion, which was confirmed by HE AO 1. The flux 
measurements of 4U 1543—624 are roughly constant 
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odicities from 50 s - 10,000 s found in the EXOSAT 
data, and no periodicitie s from 0. 1 s - 1 000 s detected in 
the SAX data. Recently, ISchultd l)2002|) presented spec- 
tral results from archival ASCA, SAX, and RXTE ob- 
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servations. These results suggest that an increase in the 
luminosity of 4U 1543—624 is accompanied by a hard- 
ening of the spectrum. In addition, an Fe-K emission 
line is see n in the hard s tate, but not detected in the 
low state. iFarinelli et"afl {2003) present an independent 
analysis of the two SAX observations which show no sig- 
nificant luminosity change b ut a spectral h a rdenin g in 
the second obser vation. Both IFarinelli et al.l l)2003|) and 
lAsai et alJ lj2000ft report Fe-K line detections for the SAX 
and ASCA observations, respectively. 

In contrast, 2S 0918—549 shows a factor of 10 X-ray 
variability, but again has no known periodicities in either 
the X-ray or optical bands down to timescales of 1 hour 
(Forman et al. 1978; Warwick et al. 1981; Chevalier 
& Ilovaisky 1987; Smale & Lochner 1992; Christian & 
Swank 1997; Schulz 1999; Jonker et al. 2001). Cheva- 
lier & Ilovaisky (1987), using the Einstein HRI position, 
identified an ultraviolet-bright optical counterpart for 
2S 0918—549 and suggested a source distance of 15 kpc 
based on the properties of other LMXBs, i.e., My=0.0 
and (B - V) =0.0. Recently, Jonker et al. (2001) de- 
tected a type I X-ray burst from 2S 0918—549 and de- 
rived an upper limit to the distance of 4.9 kpc from ra- 
dius expansion arguments. A similar analysis using Bep- 
poSAX Wide Field C amera data of 2S 0918- 549 implied 
a distance of 4.2 kpc (ICorncliss e et al.ll20 02'). 

We have identified both 2S 0918-549 and 
4U 1543-624 as being part of a class of four NS/LMXBs 
all having a simila r feature at 0.7 k eV in their low- 
resolution spectra Ipuett et alJ 1200 lj) . This feature 
had been attributed to unresolved line emission from 
Fe and O (see, e.g., iChristian. White fc Swank! 119941 
White, Kallm an fc Angelinil I1997JK However, a high- 
resolution observation of the brightest of these sources, 
4U 0614+091, with the Chandra X-Ray Observatory, 
failed to detect any emission lines, finding instead an 
unusually high Ne/O number ratio in the absorption 
along the line of sight (Paerels et al. 2001). Previously, 
we showed that the ASCA spectra of all four sources are 
well fit without a 0.7 keV emission line using a model 
that includes photoelectric absorption due to excess Ne 
along the lines of sight and presumably local to the 
sources. In this paper, we present the high-resolution 
Chandra and XMM-Newton spectra of two sources, 
2S 0918-549 and 4U 1543-624, which support our 
earlier ASCA results. 

Given the high Lx/L opt ratio of these sources and the 
excess Ne absorption, we suggested that the systems are 
ultracompact bin aries with Ne-rich degenerate donors 
ijJuett et al.M200i|) . This is similar to the ultracompact 
X-ray pulsar 4U 1626—67. We note that while the spec- 
trum of 4U 1626—67 does show O and N e emission lines 
ijAngelini et alJll995USchulz et alJl2001|) . it does not re- 
semble the sources we are considering in either low or 
high-resolution spectra. While we do compare our re- 
sults to those of 4U 1626—67, it is important to consider 
the differences between the systems when drawing con- 
clusions. 

2. OBSERVATIONS AND DATA REDUCTION 

Chandra observed 2S 0918-549 on 2000 July 19 and 
4U 1543-624 on 2000 September 12 for 30 ks each us- 
ing the High Energy Transmission Grating Spectrometer 
(HETGS) and the Advanced CCD Imaging Spectrome- 



ter (ACIS; Weisskopf et al. 2002). The HETGS carries 
two transmission gratings: the Medium Energy Grat- 
ings (MEGs) with a range of 2.5-31 A (0.4-5.0 keV) and 
the High Energy Gratings (HEGs) with a range of 1.2- 
15 A (0.8-10.0 keV). The HETGS spectra are imaged by 
ACIS, an array of six CCD detectors. The HETGS/ACIS 
combination provides both an undispersed (zeroth or- 
der) image and dispersed spectra from the gratings. The 
various orders overlap and are sorted using the intrin- 
sic energy resolution of the ACIS CCDs. The first-order 
MEG (HEG) spectrum has a spectral resolution of AA = 
0.023 A (0.012 A). 

The "level 1" event files were processed using the CIAO 
v2.2 data analysis package. For both sources, the stan- 
dard CIAO spectral reduction procedure was performed. 
During our analysis, it was found that the pipeline tool 
acis_detect_af terglow rejects 3-5% of source photons 
in grating spectra. Afterglow is the residual charge left 
from a cosmic-ray event which is released over several 
frames and can cause a line-like feature in a grating 
spectrum. The acis_detect_af terglow tool flags events 
that occur at the same chip coordinates in consecutive 
frames. For bright sources, the tool also rejects a frac- 
tion of the source events. Although only a small frac- 
tion of the total, the rejection of source photons by this 
tool is systematic and non-uniform. Since order-sorting 
of grating spectra provides efficient rejection of back- 
ground events, the afterglow detection tool is not nec- 
essary. Therefore, we reextracted the event file, retain- 
ing those events tagged by the acis_detect_af terglow 
tool. No features were found that might be attributable 
to afterglow events. 

The combined MEG + HEG first order dispersed spec- 
trum of 2S 0918—549 has an average count rate of 
3.10±0.18 cts s _1 . We examined the total count rate, as 
well as the count rates in three different energy ranges, 
to check for changes in the spectral state of 2S 0918—549. 
We found no evidence for any change of state during the 
Chandra observation. For 4U 1543—624, the first order 
MEG + HEG count rate is 17.2±0.4 cts s _1 . Similarly, 
no changes of state of 4U 1543—624 were found. To look 
for periodic modulations of the X-ray flux, we created 
lightcurves from the combined first order event files after 
barycentering and randomizing the event arrival times. 
Randomizing of the event arrival times consists of adding 
a random quantity uniformly distributed between 0-3.2 s 
in order to avoid aliasing caused by the readout time. 
We searched for modulations of the X-ray flux with fre- 
quencies between 4xl0 -5 and 5xl0 -2 Hz. We found no 
evidence for periodic modulation in cither source, with 
a 90%-confidence upper limit of 2% (2S 0918-549) and 
1% (4U 1543—624) for the fractional rms amplitude. 

For both observations, we applied the aspect offset cor- 
rection available at the Chandra website 2 . Using the 
CIAO tool celldetect, the zeroth order source position 
for 2S 0918-549 was determined: R.A.=09 h 20 m 26!95 
and Dec=-55°12'24"7, equinox J2000.0 (90% confi- 
dence error of 0'.'6). S ince the optical counterpar t 
position is not given by iChevalier fc Ilovaiskvl l)1987|) . 
we obtained optical images of the 2S 0918—549 field 
using the Magellan 6.5-m telescopes in Chile. Wc 
measured the position of the optical counterpart by 

2 See http://asc.harvard.edu/cal/ASPECT/celmon/index.html 
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matching about 20 field stars wi th the USNO-A2 
catalog of astrometric standards ijMonet et alJ 119981) . 
We obtained R.A.=09 h 20 m 26 1 !99 and Dec=-55°12'24"8 
(equinox J2000.0) with an error radius of 1"8 (90% 
confidence). The optical and X-ray positions for 
2S 0918—549 are l'/5 separated. The zeroth order im- 
age of 4U 1543—624 was so heavily piled up that no 
counts were recorded in the center of the point spread 
function. Celldetect was not able to determine an ac- 
curate source position, but the CIAO tool tgdetect, 
which runs celldetect in conjunction with grating spe- 
cific filters and sub-tools, was able to position the zeroth 
order image of 4U 1543-624: R.A.=15 h 47 m 54!69 and 
Dec=-62°34'05'.'4, equinox J2000.0. This position is Of.' 2 
from the optical counterpart position given by Bradt & 
McClintock (1983). 

For bright sourc es, pileup ca n be a problem for CCD 
detectors (see, e.g.. lDavisl2003|) . Pileup occurs when two 
or more photons are incident on the same pixel during the 
3.2 s readout time of the ACIS detector. When this hap- 
pens, the instrument reads just one event at an energy 
comparable to the sum of the original photons energies, 
thus reducing the detected count rate and modifying the 
spectral shape. The zeroth order Chandra spectra of 
both 2S 0918-549 and 4U 1543-624 were heavily af- 
fected by pileup (pileup fraction >80%) and were not 
used in this analysis. In addition, pileup can also af- 
fect dispersed spectra. We checked the dispersed spec- 
tra of both sources and found that for 2S 0918—549 no 
pileup occurred; whereas for 4U 1543—624, there were 
signs of pileup in the first order MEG spectrum. In dis- 
persed spectra, pileup can affect only a limited wave- 
length range, particularly where the effective area of the 
instrument is the highest. For 4U 1543—624, pileup oc- 
curred in the range 5-12 A (1-2.5 keV). 

For 2S 0918-549, we used the standard CIAO tools to 
create detector response files (ARFs) for the MEG and 
HEG +1 and —1 order spectra. These were combined 
when the +/— order spectra were added for the MEG 
and HEG separately. The data were binned to 0.03 A 
(0.015 A) for the MEG (HEG) to ensure good statistics. 
We also created background files for the MEG and HEG 
spectra using the standard CIAO tool. 

There are known differences in the quantum efficiencies 
(QEs) of the front-illuminated (FI) CCDs compared to 
the back-illuminated (BI) CCDs. An initial analysis us- 
ing the standard QE files revealed that, for 4U 1543—624, 
the chip differences were identifiable (especially at high 
wavelengths) and needed to be corrected. Using a correc- 
tion table 3 provided by H. L. Marshall of the HETGS in- 
strument team, we modified the QE files available in the 
standard CIAO calibration database and then used them 
to create ARFs for the spectral analysis of 4U 1543—624. 
This modification was not required for the 2S 0918—549 
data. In order to use the grating pileup kernel in ISIS, we 
created ARFs for each chip using the CIAO tool mkgarf 
and then combined them to create +1 and —1 order 
MEG and HEG ARFs using a custom tool developed 
by J. Davis of the HETGS instrument team. (This tool 
is similar to the standard CIAO tool but also correctly 
calculates the fractional exposure at each response bin. 

3 Documentation and table available at 
http: / / spa ce.mit.edu/ CXC /calib /hetgcal.html | 



This functionality has been added to the standard tools 
in CIAO v2.3.) The grating pileup kernel models the 
effect of pileup on grating spectra, similar to the pileup 
mod el available for CCD spectra in ISIS and XSPEC 
(see, !Davisll2"003l for a discussion of grating pileup mod- 
eling). In addition, the data and responses were rebinncd 
to 0.033 A (0.017 A) for the MEG (HEG) to reflect the 
size of an ACIS event detection cell (3 CCD pixels). We 
used this binning for the spectral analysis. All Chandra 
spectral analysis of 2S 0918-549 and 4U 1543- 624 was 
performed using ISIS l|Houck fc DeNicolall200^ . 

XMM observed 2S 0918-549 on 2001 May 5 for 
40 ks and 4U 1543-624 on 2001 Feb 4 for 50 ks. 
XMM carries three different instru ments, the European 
Photon Imaging C ameras (EPICs; Striider et alJ l200ll 
iTirner e l; al^Q Ol )) the Reflecti on Grating Spectrometers 
(RGS s:lden Herder et al.l200lD . and the Optical Monitor 
(OM; iMason et alJl2001^ . The OM was not used in this 
analysis. 

The EPIC instruments consist of 3 CCD cameras, 
MOS 1, MOS 2, and pn, which provide imaging, spec- 
tral and timing data. For both observations, the pn and 
MOS 1 cameras were run in timing mode which pro- 
vides high time resolution event information by sacrific- 
ing 1-dimension in positional information. During the 
2S 0918—549 observation, the MOS 2 camera was run 
in full-frame mode, while for the 4U 1543—624 obser- 
vation, it was run in small-window mode. The EPIC 
instruments provide good spectral resolution (Ai?=50- 
200 eV FWHM) over the 0.3-12.0 keV range. There are 
two RGSs onboard XMM which provide high-resolution 
spectra (A/AA=100-500 FWHM) over the 5-38 A range. 
The grating spectra are imaged onto CCD cameras sim- 
ilar to the EPIC-MOS cameras which allows for order 
sorting of the high-resolution spectra. 

The XMM data were reduced using SAS version 5.3. 
Standard filters were applied to all XMM data and times 
of high background were excluded. The EPIC-pn data 
were reduced using the procedure epchain. Source and 
background spectra were extracted using filters in the one 
spatial coordinate, DETX. A ready-made response file 
for pn timing mode data can be found on the XMM web- 
site 4 . We reduced the EPIC-MOS data using emchain. 
Since no response file was available for MOS timing mode 
observations, the MOS 1 data were not used in this anal- 
ysis. The MOS 2 data were found to have considerable 
pileup. To reduce the effect of pileup on the spectra, we 
extracted source spectra using an annulus to excise the 
core of the point spread function where pileup is preva- 
lent. The inner 20" were excluded for 2S 0918-549 and 
the inner 35" for 4U 1543—624. Background spectra were 
extracted using an annulus centered on the source but ex- 
cluding any other point sources in the field of view. Re- 
sponse files were generated for the MOS 2 spectra using 
rmf gen and arf gen. The RGS data were reduced using 
rgsproc, which produced standard first order source and 
background spectra and response files for both RGSs. We 
grouped the pn and MOS spectra to oversample the en- 
ergy resolution of the CCDs by no more than a factor 
of three. The RGS spectra were grouped to ensure that 
each bin had at least 20 counts. All XMM spectral anal- 
ysis of 2S 0918-549 and 4U 1543-624 was performed 

4 Available at http://xmm.vilspa.esa.es/ccf/epic/ 
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using XSPEC rnau 

The 2S 0918—549 observation had average count rates 
of 81.06±0.05, 13.03±0.02, and 8.928±0.016 cts s" 1 for 
the pn, MOS 2, and combined RGS, respectively. Using 
the pn data, we checked for changes in the spectral state 
by examining the total count rate, as well as the count 
rates in three different energy ranges. There was no evi- 
dence for a change in the spectral state of 2S 0918—549 
during the XMM observation. The average count rates 
during the 4U 1543-624 observation were 238.89±0.07, 
46.88±0.03, and 23.07±0.02 cts s" 1 for the pn, MOS 
2, and combined RGS, respectively. Again, no change 
in the spectral state of 4U 1543—624 was found during 
the XMM observation. While the pn timing mode al- 
lows for a time resolution of 0.03 ms, searches for pe- 
riodic and quasi-periodic modulation at high frequency 
are hampered by a 167 Hz instrumental signal (and mul- 
tiple harmonics). Due to the low count rates of these 
observations (<Cl count per readout), the low frequency 
end of the power spectrum is affected by instrumental 
noise (see, iKuster et al.l 120021 and references therein). 
To remove the noise, we calculated the average power 
at each frequency using an 11 bin moving average and 
then divided by this value. We verified that the powers 
were distributed correctly and then searched for modula- 
tion of the X-ray flux with frequencies between 2 x 10~ 5 
and 0.5 Hz. No modulation was detected with a 90%- 
confidcnce upper limit on the fractional rms amplitude 
of 0.3% (2S 0918-549) and 0.14% (4U 1543-624). 

3. SPECTRAL RESULTS 
3.1. 2S 0918- 549 

To determine the appropriate continuum model, we 
fit the Chandra HEG and MEG spectra jointly with 
both an absorbed power-law and an absorbed power- 
law + blackbody model over the wavelength ranges 1.5— 
13 A for the HEG and 1.8-25 A for the MEG. In 
all cases, we used the in terstella r abs orption model of 
iWilms. Allen, fe McCravl ftbabs: l2000l) . In addition, we 
included an edge at 43.7 A to accoun t for the instru- 
ment al C buildup on the ACIS CCDs 5 l)Plucinskv et all 
2002). We calculated that the appropriate optical 
depth for the instrumental C edge at the time of the 
2S 0918—549 observation was 0.79 from calibration data 
points provided by H. L. Marshall. This edge is included 
in all Chandra spectral fits of 2S 0918—549. The power- 
law + blackbody model produced a significantly better 
chi-squared value than the power-law alone (significance 
greater than 99.99% as calculated by an F-test), so we 
take a power-law + blackbody model as the appropriate 
continuum model for the Chandra spectra. The best-fit 
parameter values are given in Table 1. 

We showed that the lower-resolution ASCA spectra 
of 2S 0918-549 could be well fit by including photo- 
elec tric absorption w ith non-solar abundances of O and 
Ne l|Juett et al.ll2001^) . High-resolution spectroscopy can 
test the accuracy of this hypothesis by directly mea- 
suring the edge depths, and therefore column densities, 
of neutral O, Fe, and Ne. We can then compare the 
abundance ratios to those expected for the interstellar 
medium (ISM). To determine the columns of each of 



these elements, we replaced the absorption model with 
the compa rable variable abun dance absorption model 
ftbvarabs; IWilms et al.l l200Cfl and set the abundances 
of O, Fe and Ne to zero. The Njj parameter was allowed 
to vary while all other absorption model parameters were 
frozen. The Ne and O-K edges were each modeled by a 
multiplicative edge model of the form: 



M(E) 



for E <E t 



edge j 



exp [-r(E/E edge ) 3 ] for E > E cdgc , 



(1) 



with r and E edge as variable parameters. Individual edge 
models are used since the edge positions and shapes in 
the standard absorption models are not appropriate for 
high-resolution spectral analysis (see, Wilms et al. 2000; 
Juett et al. 2003, in preparation). A Gaussian line was 
included to model the interstellar absorption line due to 
atomic O at 23.5 A. The more complex Fe-L edge at 
17.5 A was modeled using a table model, created from 
the optical constant measurements by Kortright & Kim 
(2000), which allowed for both the column density of Fe 
and the edge energy to vary. The continuum was again 
described by a power law + blackbody. The resulting 
parameters for the power-law and blackbody models are 
consistent with the low-resolution fit. The results of the 
fit can be found in Table 1 and are shown in Figure 1. 

From the absorption edge depths, the absorbing col- 
umn densities of O, Fe, and Ne can be determined 
along with the equivalent hydrogen column for each el- 
ement assuming interstellar abundances. These results 
ar e summarized in Table 2. The at omic cross-sections 
of IHenke. Gullikso n. fc Dav is (1993) were used for O 
and Ne, while the Fe-L cross-section was taken from 



iKortright fc Km] ll200l. W e use the ISM abundances 
given bv IWilms et alJ l)2000|) to calculate equivalent hy- 
drogen column densities. The wavelengths of the edges 
are consistent with measurements of neutral ISM edges 
in other LMXBs (Juett et al. 2003). We find a Ne 
column density of N^ c = (9.9 ± 1.7) X 10 17 cm -2 and 
an O column density of N = (1.9 ± 0.3) x 10 18 cm -2 
(all errors quoted at 90%-confidence level). The equiva- 
lent hydrogen column implied by the Ne edge is roughly 
three times that calculated from the O edge. We find 
a Ne/O number ratio of 0.52±0.12, compared to the 
ISM ratio of 0.18. The measured Fe column is Ap = 
(IO+3) x 10 16 cm~ 2 . The atomic O absorption line has 
a wavelength of 23.52±0.04 A and an equivalent width 
(EW) of 0.017±0.010 A. 

The evidence for Ne enhancement in the absorption to 
2S 0918-549 is sim ilar to what is seen in 4U 1626-67 
ijSchulz et alJl200Hl . Since the spectrum of 4U 1626-67 
also shows emission lines from Ne and O, we per- 
formed a careful search of the Chandra spectral resid- 
uals of 2S 0918—549 to place limits on the presence of 
any spectral line features. Gaussian models with fixed 
FWHM=2000 km s _1 were fit at each point in the wave- 
length range 1.5-13 A for the HEG and 1.8-25 A for the 
MEG. From this, we can place a 3<r upper limit on the 
EW of any line feature, either emission or absorption. 
The EW limit increases with wavelength, varying from 
0.02 A at 1.5 A to 0.5 A at 23 A. These EW limits are 
3-10 times lower than the line EWs found in the Chan- 
5 more information see, ^£ a observation of 4U 1626—67. Since the lines found in 

|http://cxc.harvard.edu/cal/Links/Acis/acis/Cal-prod s/qeDeg/index.htnthe Spectrum of 4U 1626—67 were broader, FWHM of 
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Table 1. Spectral Fits 3, 
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3.3+0.5 


0.55+0.04 


23+9 


1.07 


XMM tbabsx PL 


2.95+0.02 


2.265+0.004 


6.23+0.03 






1.61 


XMM full model 


3.20+0.03 


2. 25 (fixed) 


5.48+0.06 






1.27 






4U 1543- 


-624 








Chandra tbabsx (PL+BB) 


1.58+0.11 


1 539 +u ' ua 
i.ooa_ 004 


12.6+0.6 


0.62+0.03 


61+13 


1.75 


Chandra full model 


3.31+0.014 


1.97(flxed) 


17.1+0.8 






1.32 


XMM tbabs x (PL+BB) 


3.478+0.017 


2.568+0.008 


25.12+0.11 


1.480+0.005 


3.38+0.07 


4.34 


XMM full model 


3.613+0.014 


2.7(flxed) 


16.26+0.06 


1.512(fixed) 


18.5+0.3 


1.87 



a All errors are quoted at the 90%-confidence level 

b PL=power law; BB=blackbody; tbabs=photoelectric absorption model of Wilms ct al. (2000); full model=absorption of O, Fe, and Ne 
modeled by edge models with O absorption lines where needed. The results for the edge models are given in Table 2. 
c Power-law normalization at 1 keV in units of 10 — 2 photons keV -1 cm -2 s~ 1 
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Fig. 1. — (upper panels) Chandra MEG first order spectrum of 
2S 0918—549 with best-fit power-law + blackbody model including 
the O, Fe, and Ne edge models. For clarity, the data is binned to 
0.06 A. (lower panels) Fractional residuals ([data— model] /model) 
of the MEG first order spectral fit shown above. There are no 
features detected which are consistent with narrow emission or ab- 
sorption lines. The arrows mark the Ne-K and Fe-L absorption 
edges at 14.3 and 17.5 A respectively. 



5340 and 7460 km s" 1 for the Ne X and O VIII lines fit 
by single Gaussians, we investigated how the EW limit 
changes as the model FWHM is increased. Increasing 
the FWHM of the line to 5000 (7000) km s _1 , increases 



the EW upper limits by a factor of 1.5 (2.0). 

To determine the best-fit continuum model for the 
XMM spectra of 2S 0918-549, we simultaneously fit the 
pn, MOS 2, and both RGS spectra. Since the pn tim- 
ing mode response is not well calibrated below 1 keV, we 
only fit the pn spectrum in the 1.0-12.0 keV range. The 
energy range used for the MOS 2 was 0.3-8.0 keV and for 
the two RGS spectra it was 0.35-2.0 keV. We again fit 
both an absorbed power-law and an absorbed power-law 
+ blackbody model. We also included a multiplicative 
constant to account for any instrumental normalization 
differences. During the XMM observation, the black- 
body component was not significantly detected and we 
choose an absorbed power-law as the appropriate contin- 
uum model. The best-fit parameter values for the XMM 
spectral fit are given in Table 1. 

For the high-resolution measurements of the absorp- 
tion edge depths, we fit only the 2 RGS spectra using 
a power-law as the continuum model. Since the RGS 
covers only a limited wavelength range, 6.2-37.6 A, wo 
froze the power-law index to the best-fit value of the 
continuum model fit. Again we used the variable ab- 
sorption model with O, Fe and Ne abundances set to 
zero, and included the appropriate edge models for O, 
Fe, and Ne. A Gaussian model was used to model the 
atomic O absorption line. The best-fit results are given 
in Tables 1 and 2 and shown in Figure 2. From the fits 
of the XMM RGS spectra, we find a Ne column density 
of Vno = (9.1 ± 0.6) x 10 17 cm~ 2 and an O column den- 
sity of V = (1.96 ± 0.07) x 10 18 cm" 2 , both consistent 
with the Chandra measurements. We find a Ne/O num- 
ber ratio of 0.46±0.03. The best-fit atomic O absorption 
line had a wavelength of 23.504±0.013 A and an EW of 
0.020±0.005 A. The initial best-fit Fe-L column density 
left large residuals around the edge shape suggesting that 
the value was too large. It is possible that fits of the en- 
tire RGS wavelength range are not as sensitive to the 
Fe-L complex. To determine the correct Fe column, we 
fit only the 16.5-18 A range with a simple model con- 
sisting of a power-law and the Fe-L table model. The 
power-law index was allowed to vary for this fit. From 
this fit, we found a Fe column of (8 ± 3) x 10 16 cm~ 2 , 
which was a much better visual match and consistent 
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Fig. 2. — (upper panel) XMM combined RGS first order spec- 
trum of 2S 0918—549 with best-fit power-law model including the 
O, Fe, and Ne edge models. For clarity, the data is binned to 
0.09 A. (lower panel) Fractional residuals ([data— model] /model) of 
the RGS first order spectral fit shown above. The absorption line 
at 23.5 A is the expected atomic O line. There are no other features 
detected which are consistent with narrow emission or absorption 
lines. The arrows mark the Ne-K, Fe-L and O-K absorption edges 
at 14.3, 17.5, and 23 A respectively. 



with the Chandra value. Our values for other parame- 
ters were determined with the Fe column density fixed 
to the above value. 

We also searched for line features in the XMM spec- 
trum of 2S 0918—549 using the same procedure as 
for the Chandra spectra. Gaussian models with fixed 
FWHM=2000 km s _1 were fit at each point in the wave- 
length range 6.25-37.6 A for both RGS spectra. We place 
a 3cr upper limit of 0.05 A on the EW of any line feature 
in the XMM spectrum of 2S 0918-549 in the wavelength 
range 6-18 A. At larger wavelengths, the EW limit in- 
creases to 0.6 A at 35 A. Again, these EW limits are 
much lower than found in the 4U 1626—67 spectrum. If 
the FWHM of the Gaussian is increased to 5000 km s _1 , 
then the EW limit increases by roughly a factor of 2. 
There are residuals in the fit around 14 A and 19 A 
at the 20% level (see, Figure 2). The residuals of the 
4U 1543—624 RGS spectral fit are similar (see, Figure 
5), while the Chandra fit of 2S 0918—549 does not show 
the same features (see, Figure 1). This suggests that the 
residuals are due to remaining uncertainty in the instru- 
ment calibration. 

The absorbed 2-10 keV flux of 2S 0918-549 was 
8.1 x 10~ n and 1.1 x 10~ 10 erg cm -2 s _1 for the Chandra 
and XMM observations, respectively. This is an order of 
magnitude lower than the ASCA observation in 1995. 
The Chandra and XMM flux measurements are, how- 
ever, consistent with other past observations, including 
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Fig. 3.— The 2-10 keV absorbed flux of 2S 0918-549 in units 
of 10~ 10 ergs cm -2 s — 1 over the past 30 years. The majority 
of the observations are reasonably consistent, except for the 1995 
ASCA observation, when the source seems to have been in a bright 
state. The large diamonds represent average flux measurements 
for monitoring missions, while the small squares represent pointed 
observations. 



the RXTE ASM flux measurements made contempora- 
neously (see, Figure 3). The spectral shape of the source 
has been roughly constant over the past 20 years, ex- 
cluding the addition of a blackbody component in some 
more recent data. Although there is some difference in 
the best-fit slope of the power-law in each observation 
(r = 1.73 for ASCA, 2.04 for Chandra, and 2.265 for 
XMM), this difference is not enough to account for the 
observed change in flux. Rather, the source seems to be 
brighter during the ASCA observation with larger nor- 
malizations measured for both the power-law and black- 
body components. We suggest that 2S 0918—549 was in a 
bright state during the ASCA obs ervation, similar to out- 
burst s seen in other LMXBs (e.g., Homan & van der Klisl 
2000). As seen in X-ray color-color diagrams for sources 
in outburst, the hard color varies between outburst and 
quiescent pha ses while the soft color remains roughly 
the same (see. iMuno et alJl2002j> . ASCA, Chandra and 
XMM have energy sensitivity which cover only the soft 
color energy range in standard color-color diagrams, so a 
spectral difference during an outburst may not be mea- 
surable with these instruments. 

3.2. 4U 1543-624 

We fit the four (+1 and -1; MEG and HEG) Chan- 
dra spectra of 4U 1543—624 jointly with both an ab- 
sorbed power-law and an absorbed power-law + black- 
body model in the wavelength range 1.7-25 A for the 
MEG and 1.3-13 A for the HEG. The pileup kernel in 
ISIS was used to model the pileup in the grating spec- 
tra. In addition, an edge model at 43.7 A was included 
to account for the instrumental C edge. We calculated 
that the appropriate optical depth for the instrumental 
C edge at the time of the 4U 1543—624 observation was 
0.85. This edge is included in all Chandra spectral fits 
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of 4U 1543—624. The absorbed power-law + blackbody 
model was a significantly better fit to the data and we 
take this as the appropriate continuum model. The ab- 
sorbed 2-10 keV flux of 4U 1543-624 during the Chan- 
dra observation was 7.1xl0 -10 erg cm~ 2 s _1 , which is 
consistent with past observations. 

As was done with the 2S 0918—549 observations, we 
modeled the O and Ne edges with edge models and the 
Fe-L edge w ith the table mode l based on the latest cross- 
section from lKortright fe Kiml (J2000) . A Gaussian model 
was included to fit the atomic O absorption line at 23.5 A. 
Since the most interesting features in the spectrum are 
at > 12 A, where pileup is no longer a problem, we fit the 
MEG spectra in the wavelength range 12-25 A without 
the pileup kernel to measure the edge features. Pileup 
in the grating spectra enhances the instrumental resid- 
uals found at 2 keV, inflating the chi-squared value of 
the fit. By fitting only the high wavelength end of the 
spectra, we determine a more characteristic chi-squared 
value which can be used to judge the goodness of the 
fit. The blackbody component is only prominent at low 
wavelengths, so we chose to model the continuum in this 
limited wavelength range as a power-law. We found that 
if the equivalent hydrogen column density of the absorp- 
tion model and the power-law index are allowed to vary, 
there is a degeneracy between the parameters that drives 
./Vh and T to higher values. To remove this degeneracy, 
we fixed the photon index to that found in a continuum 
fit with the full model. The results can be found in Table 
2. 

From the absorption edge depths, we find a Ne column 
density of A Nc = (12.4 ± 1.1) x 10 17 cm" 2 and an O 
column density of N = (8.5 ± 1.4) x 10 17 cm~ 2 . The 
equivalent hydrogen column implied by the Ne edge is 
more than eight times that calculated from the O edge. 
We find a Ne/O number ratio of 1.5±0.3, compared to 
the ISM ratio of 0.18. The atomic O absorption line was 
found at 23.504±0.015 A with an EW of 0.04±0.02 A. 
The measured Fe column is Ap = (3 ± 2) x 10 16 cm~ 2 . 

As can be seen in Figure 4, no emission or absorp- 
tion lines are apparent in the spectrum of 4U 1543—624. 
We performed a careful search of the Chandra spectral 
residuals to place limits on the presence of any spectral 
features in the 12-25 A wavelength range of the MEG. 
Gaussian models with fixed FWHM=2000 km s _1 were 
fit at each point. From this, we can place a 3er upper 
limit on the EW of any line feature, either emission or 
absorption. The EW limit increases with wavelength, 
varying from 0.02 A at 12 A to 1.0 A at 23 A, while 
broader line models give EW upper limits a factor of 2 
larger. The EW limits found for the Chandra spectrum 
of 4U 1543—624 are lower than the emission line EWs 
found in the spectrum of 4U 1626—67. There have been 
a number of reports of a Fe K line detection in the spec- 
trum of 4U 1543—624. There is no apparent line in the 
Fe K emission region (6.4-6.7 keV) of the Chandra spec- 
trum. To place a limit on the EW of such a line, we fit 
an absorbed power-law + blackbody + Gaussian model 
to the full spectrum with a fixed line energy of 6.4 keV. 
The upper limit on the EW of the line is 60 eV. 

The analysis of the XMM spectra of 4U 1543-624 is 
similar to the XMM analysis of 2S 0918-549. We first 
jointly fit the EPIC-pn, MOS 2, and RGS spectra with 
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Fig. 4. — (upper panel) Chandra combined MEG first order 
spectrum of 4U 1543—624 with best-fit power-law model including 
the O, Fe, and Ne edge models. For clarity, the data is binned to 
0.06 A. (lower panel) Fractional residuals ([data— model] /model) of 
the MEG first order spectral fit shown above. The absorption line 
at 23.5 A is the expected atomic O line. There are no other features 
detected which are consistent with narrow emission or absorption 
lines. The arrows mark the Ne-K, Fe-L and O-K absorption edges 
at 14.3, 17.5, and 23 A respectively. 



both an absorbed power-law and an absorbed power- 
law + blackbody. For 4U 1543—624, the power-law + 
blackbody was the best-fit continuum model. During 
the XMM observation, the absorbed 2-10 keV flux was 
4.2xl0 -10 erg cm~ 2 s _1 which at the lower end of the 
range of observed flux of the source. The best-fit pa- 
rameters of the XMM spectral fit of 4U 1543-624 has a 
steeper power-law photon index than found in the Chan- 
dra observation of the source. In addition, the blackbody 
temperature and normalization are different, although 
this difference could be due to incorrectly accounting for 
the absorption features. If the full Chandra spectrum is 
fit with a model that allows for a non-ISM Ne/O ratio, 
the blackbody temperature and normalization become 
consistent with the XMM result, which itself is consis- 
tent with previous results (see, e.g.. iSchultzl 12002'). 

We added a Gaussian feature centered at 6.4 keV 
to the continuum model and determined the best fit 
width and flux of the line. We find that a wide, 
FWHM=2.64±0.14 keV, is allowed with a best-fit EW 
of 370±40 eV. This line is comparable to that found by 
Schultz (2002) in the RXTE data of 4U 1543-624. If 
the energy of the Gaussian is allowed to vary, the Gaus- 
sian model becomes even wider, FWHM«6 keV, with a 
central energy of 4 keV. While it is interesting to note 
that both the XMM and RXTE data shows this similar 
feature, the large width and instability of the XMM fits, 
leads us to believe that the feature is not a detection of 
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Fe line emission. 

For the high-resolution measurements of the absorp- 
tion edge depths, we fit only the 2 RGS spectra using 
a power-law + blackbody as the continuum model. We 
froze both the power-law index and the blackbody tem- 
perature during the fits. Again we used the variable ab- 
sorption model with O, Fe and Ne abundances set to 
zero, and included the appropriate edge models for O, 
Fe, and Ne. A Gaussian model was used to model the 
atomic O absorption line and another was included to fit 
a line feature at 23.36 A which has been attributed to 
absorption by O in metallic compounds (i.e. oxides; e.g., 
Schul zet al] 120021) but is also consistent with the O II 
ls-2p transition (see, Juett et al. 2003). We also found 
that the O edge could be better fit by including 2 edge 
models. The best-fit O edge wavelength is dependent 
on the form of the O (atomic vs. metallic compounds) 
and the presence of multiple forms of O would produce 
multiple edges. 

The best-fit results are given in Tables 1 and 2 and 
shown in Figure 2. From the fits of the XMM RGS spec- 



17 



cm 



tra of 4U 1543-624, we find A Nc = (8.0±0.3)xl0 
and A Fc = (4.7 ± 0.6) x 10 16 citT 2 . For a single O 
edge, the best-fit wavelength is 23.12±0.04 A with a col- 
umn density ofN o = (14.8±0.2)xl0 17 cm" 2 . Fitting 
with two O edges yields wavelengths of 22.86±0.04 A 
and 23.19±0.03 A and column densities of (8.2 ± 0.3) x 
10 17 cm~ 2 and (6.5 ± 0.3) x 10 17 cm~ 2 . The lower wave- 
length is consistent with the predicted O edge for atomic 
O while the higher wavelengths is expected from O bound 
in molecules. We find a Ne/O number ratio of 0.54±0.03. 
For the atomic O absorption line, we find a best-fit wave- 
length 23.515±0.006 A and an EW of 0.0095±0.0011 A. 
The other line has a wavelength of 23.371±0.009 A and 
an EW of 0.0083±0.0014 A. Only the Fe-L column den- 
sity is consistent between the Chandra and XMM obser- 
vations. The XMM spectral analysis shows an increase 
in the O column density by a factor of 1.7, relative to 
the Chandra result, and a reduction by 60% in the Ne 
column density. While smaller than found in the Chan- 
dra spectrum, the Ne/O ratio of 0.54 is still significantly 
larger than the ISM value of 0.18. 

We also searched for line features in the XMM spec- 
trum of 4U 1543—624 using the same procedure as 
for the Chandra spectra. Gaussian models with fixed 
FWHM=2000 km s _1 were fit at each point in the wave- 
length range 6.2-35 A for both RGS spectra. We place a 
3cr upper limit of 0.02 A on the EW of any line feature in 
the XMM spectrum of 4U 1543-624 in the wavelength 
range 6-18 A. At larger wavelengths, the EW limit in- 
creases to 0.3 A at 35 A. These EW limits are much 
lower than found in the 4U 1626—67 spectrum. When 
the FWHM=5000 km s -1 , the EW limits are a factor of 
2 greater. As mentioned earlier, there are definite resid- 
uals in the spectral fit of 4U 1543—624 (see, Figure 5). 
We attribute these residuals to instrumental calibration 
uncertainty. 

3.3. Residual Features 

We now address the residuals found in our spectral 
analysis. For the XMM spectra of both sources, there 
are broad residual features at the 10-20% level. These 
residuals are consistent with the uncertainty in the cali- 
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Fig. 5. — (upper panel) XMM combined RGS first order spectrum 
of 4U 1543—624 with best-fit power-law + blackbody model includ- 
ing the O, Fe, and Ne edge models. For clarity, the data is binned 
to 0.09 A. (lower panel) Fractional residuals ([data— model] /model) 
of the RGS first order spectral fit shown above. The absorption line 
at 23.5 A is the expected atomic O line while the line at 23.37 A 
is consistent with O absorption in molecular compounds. There 
are no other features detected which are consistent with narrow 
emission or absorption lines. The arrows mark the Ne-K, Fe-L and 
O-K absorption edges at 14.3, 17.5, and 23 A respectively. 



bration of the RGS 6 . All spectra show some level of struc- 
ture around the neutral Ne-K edge at 14.3 A. For both 
instruments, the RGS and HETG, the interpretation of 
this feature is complicated by instrumental features at 
the same wavelength: the RGS has lost one of the CCDs 
that covered the 10.5-14.0 A wavelength range reducing 
the effective area in this range, while the HETG nominal 
pointing places a chip gap in the —1 MEG order at the 
Ne-K edge. 

While we do not claim a significant detection of any 
feature, the similar structure seen in all spectra suggest 
that the simple interpretation of a neutral Ne-K edge 
may not be entirely correct. Below we detail some possi- 
ble astrophysical origins, remembering that instrumental 
effects are likely also important. 

One explanation is that the feature is an emission line 
or a radiative recombination continuum (RRC) emission 
feature. At 14.0-14.3 A, emission from a number of dif- 
ferent ionization states of Fe is possible, but we would 
expect to see multiple features (e.g., as in the spectrum 
of Cyg X-l; Schulz ct al. 2002). The O VIII RRC feature 
is also consistent with a 14.2 A position. Again though, 
the lack of other O VIII emission features, in particular 
the O VIII emission line at 18.97 A, makes this possibility 

6 See document XMM-SOC-CAL-TN-0030 available at 
|http: / /xmm. vilspa.esa.es/external/xmm_sw_cal/ calib / index, shtml 
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Table 2. Photoelectric Absorption Edge Measurements 3, 



Edge 


A Optical depth 
(A) ' r 


N z 

(10 17 cm- 2 ) 


iV H 
(10 21 cm" 2 ) 


2S 0918-549: Chandra result 


O K 


23.0±0.3 1.08±0.18 


19±3 


3.9±0.6 


Fe L III 


17.49±0.02 b 0.7±0.4 


1.0±0.6 


3.7±2.2 


Ne K 


14.24+_£ 02 0.36±0.06 


9.9±1.7 


11.4±1.9 




2S 0918-5A9-.XMM 


result 




K 


22.99±0.08 1.11±0.04 


19.6±0.7 


4.00±0.14 


Fe L III 

Ne K 


17.50±0.03 b 0.6±0.2 
14.40±0.03 0.33±0.02 


0.08±0.03 
9.1±0.6 


3.0±1.1 
10.4±0.6 


4U 1543-624: Chandra result 


O K 


22.72+Jj;j£| 0.48±0.08 


8.5±1.4 


1.7±0.3 


Fe L III 

Ne K 


17.49±0.04 b 0.21±0.14 
14.330±0.017 0.45±0.04 


0.3±0.2 
12.4±1.1 


1.1±0.7 
14.2±1.3 




4U 1543-624:AA/A/ 


result 




O K 
O K 


22.86±0.04 0.64±0.016 
23.19±0.03 0.368±0.016 


8.2±0.3 
6.5±0.3 


1.67±0.06 
1.32±0.06 


Fe L III 


17.50±0.02 b 0.33±0.04 


0.47±0.06 


1.7±0.2 


Ne K 


14.286j:o°5 2 0.291±0.010 


8.0±0.3 


9.2±0.3 



a All errors are quoted at the 90%-confidence level. Ato mic cross-sections t aken from Henke et al. (1993) for O and Ne, and Kortright & 
Kim (2000) for Fe-L. Interstellar abundances taken from IWilms et all 1200011 , 
b We define the absorption edge energy for Fe as the energy at which the edge structure reaches its minimum value. 



less probable. 

The O VIII absorption edge is also located at 14.2 A 
and the Chandra spectrum of 4U 1543—624 (see, Figure 
4) seems to suggest that 2 edges, one at Ne-K and another 
at O VIII may be a plausible solution. Finally, the fea- 
ture might be related to the Ne absorption edge. Unlike 
other elements, Ne cannot be in molecular forms which 
might cause structure at the K shell absorption edge. We 
can not rule out however unanticipated structure in the 
edge shape, in particular if the Ne-rich absorbing mate- 
rial has some velocity, we might expect to find a velocity 
shifted edge structure. The interpretation of this feature 
is hampered by the lack of high-resolution atomic data 
on X-ray absorption edges. As we acquire more data on 
X-ray absorption edges, it will become possible to deter- 
mine if this structure is "normal" or instead suggests the 
presence of more interesting phenomena. 

4. DISCUSSION 

We have shown that the ChandrafHETGS and 
XMM/RGS spectra of 2S 0918-549 and 4U 1543-624 
are well fit by models that allow for absorption columns 
of neutral Ne and O with abundance ratios significantly 
different from the expected ISM ratio. The Chandra 
spectrum of 2S 0918—549 has a best fit Ne/O number 
ratio of 0.52±0.12, or roughly 3x the ISM value of 0.18. 
The XMM fit of 2S 0918-549 gives a consistent Ne/O 
number ratio of 0.46±0.03. The best-fit model for the 
Chandra spectrum of 4U 1543—624 has a Ne/O number 
ratio of 1.5±0.3, while the best-fit XMM model has a 
Ne/O number ratio of 0.54±0.03. The unusual abun- 
dance ratios as well as the variability seen in the obser- 
vations of 4U 1543—624 lead us to conclude that there is 
absorption local to these binaries and that the material 
is Ne enriched. 

From the L x /L opt ratios of 2S 0918-549 and 
4U 1543—624, we expect these binaries to have or- 
bital periods <60 min based on the empirical relation- 
ship determined by van Paradijs & McClintock (1994). 



We searched the Chandra lightcurves for orbital mod- 
ulation, but found no modulations larger than 1-2%. 
XMM provided a more sensitive search for orbital mod- 
ulations, but we found no signatures of orbital mod- 
ulation at fractional rms upper limits of 0.14-0.3%. 
Confirmation of such short orbital periods would place 
2S 0918-549 and 4U 1543-624 in the class of ultracom- 
pact LMXBs (P O rb^80 min). Ultracompact binaries re- 
quire H-depleted or degenerate dwarf companions (Joss 
et al. 1978; Nelson et al. 1986). Such companions would 
be expected to have non-standard abundances compared 
to ISM values. The Chandra spectrum of 4U 1626—67 re- 
vealed absorption edges of C, O, and Ne which are 5 times 
larger than would be predicted given the hydrogen col- 
umn density measured in the UV. For 4U 1626—67, the 
measured local abundance ratios, if the excess material is 
assumed to have originated around the binary, are con- 
sistent with the expected abundances in the chemically 
fractionated core of a C-0 or O-Ne-Mg WD (Schulz et al. 
2001). The evidence for Ne absorption in 2S 0918-549 
and 4U 1543—624 hints at a similarity between these 
sources and 4U 1626—67. In addition, there is evidence 
that the optical spectra of this group of sources is also 
similar to 4U 1626—67 with no H or He lines detected, 
but with a C III/N III emission line near 4640 A (Wang & 
Chakrabarty 2003, in preparation). Based on the similar- 
ity between these sources and 4U 1626—67, we previously 
attributed the excess Ne absorption in 2S 0918—549 and 
4U 1543—624 to material local to the sources, and sug- 
gested that these sy stems contained a Ne-rich degenerate 
donor ijJuett et al.ll2001fl . 

The one assumption we have made in this analysis is 
that the absorption is from neutral material only. If there 
is a sizable contribution to the absorption from mate- 
rial local to the binary, this assumption is most likely 
not valid, since we would expect ionization of the local 
material from the central source. Accounting for ioniza- 
tion of the material may explain the large Ne/O ratios 
observed. We find Ne/O ratios larger than the 0.22 in- 
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ferred for the lo cal absorption in t he Chandra spectrum 
of 4U 1626-67 l|Schulz et al.ll200l . Ne/O ratios for lo- 
cal absorption in our source will be even higher than the 
measured total line-of-sight ratio once ISM contributions 
are removed, since the O column from the ISM should 
be larger than the Ne column. Unfortunately, we can 
not do this calculation, as was possible for 4U 1626—67, 
since we do not know the expected ISM contribution for 
2S 0918-549 or 4U 1543-624. If the local material is 
affected by ionization, the effect will be different for each 
element. O will become ionized before Ne, leading to an 
enhanced Ne/O ratio as measured by the neutral edges. 
While it is tempting to assume that ionization is the only 
cause for the unusual abundances, we point out that if 
the local material was of standard abundances, we would 
expect local absorption from higher Z elements, like Mg 
and Si, which are not seen in the spectra of 2S 0918—549 
or 4U 1543—624. This leads us to conclude that material 
must have some enhancement of Ne to show such strong 
absorption. 

One of the most interesting results is the difference in 
the Ne and O column densities in the two observations of 
4U 1543-624. The XMM results show a increase in the O 
column density with a decrease in the Ne column density. 
In addition, the XMM spectrum is both softer and has a 
lower luminosity in the 2-10 keV band than the Chandra 
spectrum, which suggests that the source was not in the 
same state in the Chandra and XMM observations. We 
suggest that the difference in the 4U 1543—624 results 
provides support for the local nature of the absorption. 
Since the high-resolution results of 2S 0918—549 are con- 
sistent not only in the measured Ne/O number ratio but 
also in the spectral model in general, we are confident 
that instrumental differences do not have a significant 
effect on the spectral results. It is interesting to note 
the the Ne/O ratios found in the Chandra and XMM ob- 
servations of 2S 0918—549 are smaller than the i nferred 
value from the ASCA result of Duett et"aTI l)2001(l . while 
at the same time the flux of 2S 0918—549 has decreased 
by a factor of 10 since the ASCA observation. It is pos- 
sible that the state changes alter the ionization structure 
of the local material, which in turn changes the absorb- 
ing columns of neutral Ne and O. It will be difficult to 
separate the effect of unusual abundances from ionization 
changes without more data on these sources. We propose 
that multiple observations of 4U 1543—624 using simul- 
taneous high-resolution and broadband observations will 
probe the connection between spectral state differences 
and the column density variations. It is important to 
note that ionization effects and the ISM contribution 
must be understood before it is possible to determine 
the intrinsic Ne/O abundance ratio of the local mate- 
rial, which can then be used to place constraints on the 
composition of the companion. 

Although we attribute the excess Ne absorption in 
2S 0918-549 and 4U 1543-624 to local material in both 
of the binaries, we also consider the possibility that the 
excess Ne is due to enhancements of the ISM along the 
line of sight. The best measurement of absorption to- 
ward an X-ray binary, the Chandra spectral analysis of 
Cyg X-l, shows columns of O and Ne that are consis- 
tent with standard ISM abundances (Schulz et al. 2002). 
To quantify the variations of the ISM abundances, we 
have undertaken a study of the ISM using column density 



measurements from the spectra of X-ray binaries (Juett 
et al. 2003). Initial results show that other LMXBs do 
not show the Ne/O ratios measured for these sources. At 
most, Ne/O number ratios are only twice the expected 
ISM values. This is still significantly smaller than we find 
in 2S 0918-549 and 4U 1543-624, but does suggest that 
there are some variations between lines-of-sight. We still 
favor the interpretation that there is absorbing material 
local to the binaries, and that this is strengthened by the 
observation-to-observation variations in the Ne/O num- 
ber ratio of 4U 1543—624 which could not be due to the 
ISM. 

Analysis of the Chandra and XMM spectra of 
2S 0918-549 and 4U 1543-624 found no Ne or O lines 
like those seen in 4U 1626—67. The absence of strong 
lines in the spectra of 2S 0918-549 and 4U 1543-624 
demonstrates that their circumbinary environments are 
different compared to 4U 1626—67. One possible reason 
for this difference may be that 4U 1626—67 is a pul- 
sar with B=Z x 10 12 G (Orlandini et al. 1998), while 
2S 0918-549 and 4U 1543-624 are likely to have weak 
magnetic fields (~10 8 G). This undoubtedly results in 
a substantial difference in the accretion flow geometry. 
There is weak evidence for a line-like feature at 6.4 keV 
in the XMM spectrum of 4U 1543-624. This feature is 
very broad and could possibly be due to a mis-modeling 
of the continuum emission rather than a Fe line. In addi- 
tion, there is evidence for structure around the Ne edge 
in the spectra of both sources. It is likely that this is due 
to instrumental effects but could also be due to structure 
in either neutral Ne edges or possibly due to confusion 
with ionized O edges. 

Recent RXTE and BeppoSAX observations of 
2S 0918—549 found the first thermonuclear X-ray bursts 
from this source ijJonker et alJ 120011 iCornelisse et alJ 
2002). Both bursts were short, with durations <100 s. 
These short bursts suggest that H and/or He is undergo- 
ing unstable burning on the surface of the NS. If the com- 
panion has a large C abundance, we might expect to see a 
much longer burst («1 hr), like the "superbursts" seen in 
some LMXBs, which have been attributed to thermonu- 
clear burning of C on neutron star surfaces (Cornelisse 
et al. 2000; Cumming & Bildsten 2001; Strohmayer & 
Brown 2001). The more ordinary properties of the X-ray 
bursts from 2S 0918—549 may indicate that the donor is 
not a C-0 dwarf as suggested for 4U 1626—67. In this 
case, another mechanism for Ne enhancement is needed. 
On the other hand, it may still be possible for a car- 
bon accretor to show short X-ray bursts. One possibility 
is that while the companion is H-deficient, there is still 
a non-negligible H fraction that is accreted by the NS 
and is then respons ible for the Type I X-ray bursts (e.g., 
iNelson et alJ fl986') . Alternatively, the heavy elements 
(C, O, and Ne) may undergo spallation during accre- 
tion, leaving He and H nuclei which could then undergo 
unstable thermonuclear burning in the usual way (Bild- 
sten, Salpeter, & Wasserman 1992). Spallation reactions 
would produce 7-ray emission lines at 4.4 and 6.1 MeV. 
Unfortunately, the strength of these lines as calculated 
by Bildsten et al. (1992) is below current observational 
detection limits. For a 1 Ms observation with Integral, 
the detection sensitivities are ^1000 times higher than 
the most optimistic line flux estimates. 
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